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^ ■ Abstract 

j3 ; 

1-^ , It is shown that an (anti-)self-dual homogeneous vacuum gluon field appears 

ly^ ' in a natural way within the problem of calculation of the QCD partition function 

CN . in the form of Euclidean functional integral with periodic boundary conditions. 

I , There is no violation of cluster property within this formulation, nor are parity, 

^ I color and rotational symmetries broken explicitly. The massless limit of the 

' product of the quark masses and condensates, mf{ijjf'ijjf), is calculated to all 

. loop orders. This quantity does not vanish and is proportional to the gluon 

I condensate appearing due to the nonzero strength of the vacuum gluon field. 

' We conclude that the gluon condensate can be considered as an order parameter 

. both for confinement and chiral symmetry breaking. 

1 Introduction 

(D : 

' The physical picture of nonperturbative QCD vacuum realised with the (anti-)self-dual 

>• ■ homogeneous gluon field has become prominent since the early eighties, when Leutwyler 

^ . demonstrated the stability of this gluon configuration against local quantum fluctua- 

tions and noticed, that this field could be related to the problems of confinement and 
chiral symmetry breaking [|l], |[. Elizalde and Soto, and many other authors have ob- 
tained strong evidence that this field could be a true minimum of the QCD effective 
potential(see [|], |[ and references therein). Manifestations of this gluon configuration 
in the spectrum and weak decays of light mesons, their excited states, heavy quarkonia 
and heavy-light mesons were studied in recent papers 0, p. The vacuum field un- 
der consideration produces several qualitative regimes for masses and decay constants 
which are completely consistent with experimental data. Namely, the masses of light 
pseudoscalar and vector mesons are strongly split, orbital and radial excitations of 
light mesons show Regge behaviour, the mass of heavy quarkonium tends to be equal 
to sum of the masses of quarks, the heavy-light meson mass approaches the mass of 
the heavy quarks, and the weak decay constant for pseudoscalar heavy-light mesons 
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has asymptotic behaviour 1 / ^vtlq . Moreover, scalar and axial mesons are absent in 
the spectrum as simple qq states, but appear in the super-fine structure of orbital ex- 
citations of vector mesons. Quantitatively, the masses and decay constants of mesons 
from all different regions of the spectrum are described within ten percent inaccuracy. 
These different phenomena are displayed with the minimal set of parameters: gauge 
coupling constant, strength of the vacuum field and the quark masses. It looks as if 
the field under consideration produces both confinement and chiral (fiavour) symmetry 
breaking 0. However, there are three essential gaps that hinder justification of this 
physical picture. Regular formulation of the problem about QCD ground state, real- 
ized by the (anti)-self-dual homogeneous field, is missed. There is no proof that this 
field minimizes the QCD effective potential. Most of the results concerning the field 
under consideration are obtained within the one-loop approximation. In this paper, we 
attempt to fill in the first and third gaps. 

We construct a representation for the Euclidean QCD generating functional which 
includes the vacuum field under consideration in a self-consistent manner. Using this 
representation, we investigate the massless limit of the renormalization group invariant 
quantity m{%ljip)B which is a product of the quark mass and quark condensate in the 
presence of vacuum (anti-)self-dual homogeneous gluon field. An important result of 
this work is the formula 




where B is the strength of the vacuum field (gauge coupling constant is included into 
B). Equation (|l]) is valid to all loop orders. 

2 Generating functional 

First of all, we need to explain what is hidden behind the symbol (. . In other 
words, what is the formal statement of the problem about QCD ground state, within 
which this field appears in a natural and self-consistent way? 

The usual statement of the problem about vacuum (phase) structure of quantum 
field systems is based on the analogy between the functional integrals in Euclidean 
QFT and the partition function of quantum statistical systems in the infinite volume 
(thermodynamic) limit. Therefore, we need an appropriate representation for the QCD 
partition function in the infinite volume limit. The most subtle point is a choice of 
boundary conditions for the functional space of integration. The standard way is to 
introduce a large space-time box, to impose periodic boundary conditions on the fields 
in the box and then to study the infinite volume limit. We follow just this prescription. 

It should be noted, that condition for the fields to vanish at infinity, which is 
normal for QFT in the perturbative regime, is not appropriate, since the translation 
invariant fields are excluded ad hoc. There is no chance to get insight into the critical 
phenomena of long range correlations. The instanton-like formulation of the problem 
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- to calculate a transition amplitude from the field configuration A given at Euclidean 
time Ti to another configuration A' at time T2 - is not suitable either. In this 
case, the homogeneous field comes through the boundary conditions, which results in 
a hard violation of the cluster property and explicit breakdown of rotational and color 
symmetries and parity. 

Let us start with pure gluodynamics. A naive representation for partition function 
looks like 

DAexpf^J^d^xCYuiA)^ (2) 

where V is a large Euclidean volume, L and l3~^ = T are the space box size and the 
temperature. The functional space J-'l,/3 contains gauge fields A^ satisfying periodic 
boundary condition. Notice, that translation invariant fields, in particular the (anti- 
) self-dual homogeneous field 



Bp{x) = -nB^^x^, = tgsin^ + tgcos^, (3) 



B^u = i-B^i/, B^pBpy = —S^uB'^, B^ = const, 
belong to J^l/s- In case (|^), an arbitrary translation 

B^{x + = Bp{x) +d^uj{x,0, uj = x,B,{^), (4) 

can be compensated by a suitable gauge transformation. 

Field configuration (|^) is not a dynamical variable in the sense, that its equation 
of motion does not contain any derivatives, but is just a constraint. This field must be 
integrated out if one looks for an integral representation for partition function which 
corresponds to an actual ground state of the system. However, this integration should 
be based on resolving the constraint which takes into account all quantum corrections 
coming from the dynamical modes of the gauge fields. The quantum constraint can 
have solutions that are neither visible at the classical level nor within the perturbation 
theory. At the same time, these nontrivial solutions for the constant fields (condensates) 
govern critical phenomena in systems with the infinite number of degrees of freedom. 
One can easily illustrate this statement by the phase transitions in the models with 0^ 
and Yukawa interactions (e.g., see ^ and references therein). 

The integral over the homogeneous field can be separated in (^ with a simultaneous 
fixing of a gauge of dynamical fields by means of the Faddev-Popov trick: 

00 

1 = $[^] J DA J Duj J dB J dae <5 - - B'^] 6 [V''^{B)A''^ 
f Ss 

where the space T does not contain the nondynamical mode (|^), i? is the field strength, 
V^(-B) = — iBp[x) denotes a covariant derivative in the adjoint representation. The 
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coupling constant g is included into the field S^. The measure das is defined as 



2tt 7t 2n 



das = j^Y. I I dO sine I dC = l. (5) 







Definition (|^) corresponds to integration over the spatial (spherical) angles (<y9, 6) of 
the field @ and angle ^ which defines its orientation in color space (in the diagonal 
representation of matrix n in (^)). The sign '±' corresponds to summation of the self- 
and anti-self-dual configurations. The final representation for Z is then 



oo 

Z = Vim RaN J dae J dB J DA Afp[B, A]5 [V{B)A] 

Ss ^ 

xexpy^d^xCYM{A + B)Y (6) 

where Afp[-B,v4] is the Faddeev- Popov determinant for the background gauge con- 
dition V{B)A = 0. An appropriate regularization R\ of ultraviolet divergences and 
renormalization prescription are implied in Eq. (^). By definition, an integral over the 
fields A gives rise to an effective potential of the field B^. The background field does 
not affect general renormalizability of the theory and we rewrite Z in the form 

oo 

Z = N' J daej dBexp{-VU,s[B^;gifi),fi,/3]}, (7) 

where n is the renormalization point. As has been mentioned, the background field 
in Eq. (^ includes the coupling constant B = gB. In the background gauge, the 
composition gB is RG-invariant [^J. Furthermore, the effective potential is invariant 
under gauge and parity transformations, and space rotations (this also follows from 
the general background field method). Thus, we arrive at the expression 

oo 

Z = N' J dBexp{-VU,s[B^;g{^), 1^,(3]}. 



If the effective potential has a minimum at nonzero field strength B = B{g{fi), fi, P), 
then, in the infinite volume limit, the saddle-point method gives 

Z = exp{-VF[g{fi),fi,f3]}, 

F = U,s[B\g{fi),fi,P);g{fi),fi,P] < 0. 

The free energy density F is RG-invariant. For zero temperature, B is nothing other 
than the RG-invariant combination of the running coupling constant g{fi) and the 
renormalization point /i, hence: 

lim B' = A'qcb, limF = -CMjcD, (8) 



4 



where and Cp are positive numbers, and 

dg 



^QCD = exp I j 



These equations hnk the strength of the vacuum field with the "fundamental scale" 
Aqcd (see also Billl). 

Now we can represent the partition function Z in the form of functional integral 
over the fields A, that does not contain the translation invariant mode: 

Z = \im RaN J das J DA Afp[B, A]6 [V{B)A] exp d^xCymiA + . (9) 

Equation (P) gives the representation that we are looking for. It is based on the 
strong but single assumption that the (anti-) self-dual field corresponds to the global 
minimum of the QCD effective action. Lattice calculation of the effective potential 
for different translation invariant gluon fields seems to be the most direct way to 
verify this assumption. However, the (anti-) self-dual field is a particularly interesting 
configuration due to other reasons to be discussed below (see also [|l], H). 

On the basis of representation (|^), the QCD generating functional for correlation 
functions in the infinite volume and at zero temperature has to be defined as 

Zs[J,V,v]= lim^RAMB J das j D^ia^A.B) j J]^^/^^/ 



exp d x^jjf{x) iV — rrif + gA Tpf{x) + i J d x JA + rjfijjf + ipfrjf 

DfiA{A, B) = DAAfp [B, A]6 [V{B)A] exp | J d^xCyulA + | , (10) 

i = 7^A^, V = 7mV^, V^ = d^-iB^. 

The constant Mb provides the standard normalization Z0[O,O, 0] = 1. The functional 
space JF contains the gauge fields vanishing at infinity. The change of boundary con- 
ditions (vanishing fields instead of periodic ones) is unimportant for physics, since 
quantum fiuctuations A does not contain translation invariant modes. We have also 
added massive quarks. The fermionic functional integral spans the Grassmann algebra 
Q of square integrable fields. To be more precise, we will define this integral via a de- 
composition of the fields %() and ip over the eigenmodes ipn of the Dirac operator in the 
presence of vacuum gluon field B (an anti-hermitian representation for the 7-matrices 
in Euclidean space is used) 

- iS/i)n{x) = iXn^n{x). (11) 

As a matter of fact, this definition of the fermionic integral implies, that the ground 
state of the system is governed by the vacuum field B, and the interaction ipAip of 
quarks with the quantum gauge field A has to be treated as perturbation. Now, let us 
seek insight into the properties of representation (p^O]). 
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3 Parametrization, cluster property, symmetries 



The generating functional ([ToD contains the intrinsic dimensionful quantity B (see also 
(^), which provides the natural reference scale for running quark masses fhflfi) and 
gauge coupling constant as(/i). Therefore, the strength of the vacuum field B, the quark 
masses and coupling constant at the scale ^ = ^JB can be considered as the physical 
(intrinsic) parameters of QCD in the representation ([l0|) . Values of the parameters can 
be extracted from the analysis of hadron spectrum (e.g., see ||^). 

Correlation functions for the local or nonlocal operators Oj[A,ijj,ijj] defined in the 
standard way {A E ^ \) 



' 6 _6_ _6_' 
i6J^ idfj'' iSf] 



...Or, 



6 6 5' 
i6J^ i6f]^ i6rj 



(12) 



J=ri=ri=Q 



ensure the cluster property. Due to the integration over the angular variables 
and summation of the self- and anti-self-dual configurations, the parity, rotational and 
color symmetries are not broken explicitly. The correlators depend only on B^ and 
have normal transformation properties. At the same time, violation of the symmetries 
is seen in the integrand of Eq. (|10|), which is an indication of spontaneous breaking of 
the above-mentioned symmetries. Thus, we meet very unusual mechanism of SSB due 
to the condensation of the vector bosons. The order parameter for the nonperturbative 
phase is obvious. This is the lowest nonvanishing gluon correlator (gluon condensate), 
defined as 



) = 4i3^ + (pert.corr.) 



(13) 



Here A = {A + B) E J^, and (. . .) denotes an averaging by means of Eq. (^). However, 
any of the correlators 



{0,[A] ■ - ■ On[A]) = I dasO,[B] ■ ... ■ 0^[B] + 
contains a constant part and can be taken as order parameter. 



(14) 



4 Chiral symmetries 

Spontaneous violation of parity should influence the chiral symmetries Ua{^) and 
SUl{Nf) X SUr{Np). Due to summing the self- and anti-self-dual configurations in 
(p!OD, the vacuum expectation value of a pseudo-tensor operator is identically equal to 
zero. In particular, ?7a(1) symmetry is not broken in the sense that 

An explicit violation of parity as in the instanton ^-vacuum is needed. 
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To study the flavour chiral symmetry SUl{Np) x SUfi(Np) let us consider the 
massless limit of composition of the quark masses and quark condensates: 



Np 



/=1 



Np 



6 



6 



i6rif{x) i6rif{x 



(15) 



r]=r)=J=0 



The nontrivial point in the calculation of this quantity consists in the following. In the 
massless limit nif 0, the divergent contributions 0{mj''), with k being some positive 
integer, can appear potentially at any loop order. This means, that the perturbation 
decomposition can fail in the massless limit. In this case, the divergent terms have to 
be summed to all loop orders. This infrared problem comes /^from the zero modes ipo of 
Eq. ([TT| ) with Ao = existing due to an (anti-)self-duality of the vacuum field. It should 
be noted, that this problem arises both for the homogeneous and instanton fields. 
However, the crucial difference between the homogeneous field and the instanton 6- 
vacuum consists in the normalization of the generating functional. Unlike the instanton 
case (e.g., see [|lT], |l^), the normalization constant Af^ in Eq. ( |T0|) corresponds to 
the nonperturbative vacuum and contains the contribution of the fermion zero modes. 
Therefore, in the massless limit, no problem arises with the fermion determinant coming 
from the integral in (|10|); it is cancelled by the normalization constant A/g. 

Now we will show that a singularity 1/m exists in the lowest one-loop diagram for 
the quark condensate but does not appear in higher orders. This allows one to calculate 
the massless limit of Eq. (|15D explicitly and to prove relation (||). The most direct way 
consists in the following. 

First of all, notice that Eq. (|15D can be rewritten in the equivalent form 



Np _ d 



/=i 



/=i 



dm 



m] 



f 



Ztsim) = ^im RAAre(^) j da^ j DfiA{A,B) j H^^/^^Z 



(16) 
(17) 



exp <^ / d xipf{x) iV — mj + gA ^pf{x 



where the normalization constant N'sif^') is taken so that 2b (/i) = 1 for some yU 7^ 0. 
The LHS of Eq. (pjj) does not depend on Af^ifi), but this normalization provides us 
with an appropriately defined integral under the derivative. 

Consider for a moment the one-flavour case. An extension to Np > 1 is straight- 
forward. Formal integration over the quark field in the partition function gives 

ZB{m) = lim Rf^Afeifi) / das / DfiA{A,B)det \-iV + m~ gA 

A— >oo J J L 

Our definition of the fermion integral via the eigenmodes of Dirac operator means the 
determinant in (|I^) and its derivative in (^) are defined as 



det 



-iV + m — gA 



det 



-iV + m 



det 



gAS 



(19) 
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— det 

dm 



-iV + m — gA 



det 



-iV + m det 1 - gAS 

fr5' + -^frln(l-^i5) 
dm 



(20) 



where the trace Tr includes the space-time integration, and the quark propagator 
S{x, y) satisfies the equation 



(iV^-m^ S{x,y) = -6{x - y). 



(21) 



The term TtS in (^0]) is the lowest order contribution to the quark condensate. Higher 
perturbation corrections come from the quark loops contained in the logarithmic term 
in (^). The decisive point is that these quark loops are regular in the massless limit, 
while the lowest term is singular: 



lim IVln(l - gAS) ~ 1 + 0(m), lim TrS 

m— >0 m— >0 



1 



+ 0(1). 

m 



(22) 



It is notable, that, in another context, a regularity of the simplest two-gluon loop was 



demonstrated by Flory [13 



Using the standard representation for Green's function in terms of the matrix ele- 
ments of the projection operators Vn 



S{x,y) = 



n=0 



m + iXn ' 



one can separate the contribution of the zero eigenmodes and normal modes to the 
propagator 



S{x,y) = Soix,y) + S'{x,y), 
Soix,y) = Vo{x,y)/m, 

S'{x, y) =zV, A(x, y)P± + A(x, y) iVy + 0(m), 
^=d -iB, ^=d+iB. 

Here Vq is the projector onto the zero mode subspace 

d'^zVoix, z)Vo{z, y) = Vo{x, y), 



(23) 
(24) 



1^ 



■f{x,y)Pi 



f{x, y) = exp |-^A/n2i3(x - yY + ^-nx^B^^y^^ 



(25) 
(26) 



A{x,y) = f {x, y) / 4:TT^{x — yY is the scalar massless propagator in the background field 
(H), n is a diagonal matrix (see Eq. (H)), and P± = (1 ± 75)/2. 
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Representation (^Sj) is defined by the general square integrable solution 



V'o(a;,xo) = — — iV^uf{x,xo) (27) 

of the equation ( pTj) with Aq = 0. Details of calculation of ipQ and Vq can be found in 
Appendix. The space-time point Xq describes a position of the fermion "pseudopar- 
ticle" ipQ^x.Xo). We see that the zero eigenvalue is infinitely degenerate which is a 
manifestation of the above-mentioned invariance of the vacuum field under transla- 
tions and simultaneous gauge transformations (see (P). This feature, as well as the 
functional form of fermion zero mode (|2^), (0), is very similar to the properties of 
Leutwyler's chromons |^ which are gluon zero modes in the same background field. 

The spinor u in Eq. (^) is an eigenvector of the 75-matrix 75'?/'o = iV^o? which is 
the well-known [jTl], |14|, [1^ property of zero modes to be right-handed in a self-dual 
field and left-handed in an anti-self-dual field. As a result, the projector 



contains the chiral projection matrix P± (see (psf)). 

Representation (p^ ) for the normal mode propagator was obtained by Brown et 
al for an arbitrary ( ant i-) self- dual background field (see also ||T5|). 



Now we can return to Eq. ( ppj ) and represent the logarithmic term in the form 

(28) 



Trln(l - gAS) = J2 ^T^Tr A{So + S') 
k=i ^ 



Then one makes use of Eqs. (^) and (pSf) to notice that 

So{x,yh,A,{y)So{y,z) = (29) 

due to the projectors P± in 5*0. Therefore, all the terms in Eq. (128|), which contain 
the block Sojf^Af^So, vanish. Furthermore, any two propagators 5*0 ~ P± in the rest 
of terms of Eq. (ESf) are separated by an odd number of vertices jA and propagators 



5" ~ (7 + 0{m)). Hence, the terms in (^81) with nonzero trace of 7-matrices contain 



at least one quark mass m in the numerator for each m in the denominator. In other 
words, there is always an odd number of 7-matrices between two chiral projectors P± 
in the singular terms, and the quark loops are finite in the limit m — 0, as is pointed 
out in Eq. (^). 

Finally, taking into account equations (|16|), (|18|), ( pO]) and (p^)-(|26|), we arrive at 
lim m((ipip)i3 = — lim m — TrS* = —V^^ \ d^x^YiV^ix^x) = (30) 

m^O m^O V Jv TT 

for one fiavour. For several fiavours one gets formula ([l|). Thus, the gluon condensate 
(ITSD can be considered as an order parameter for the fiavour chiral symmetry breaking. 
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The nonzero massless limit of m{ipip)i3 indicates a non-Goldstone mechanism of 
symmetry breaking. From our point of view, breakdown of the chiral symmetry ap- 
pears here as a secondary effect of spontaneous violation of parity, which is a discrete 
symmetry. Since zero modes ( P7D are left-handed in the anti-self-dual field and right- 
handed in the self-dual field, hence in both terms (±) of generating functional (0) the 
chiral group is reduced to one of the flavour subgroups 

SUl{Nf) X SUr{Nf) SUl{Nf) (or SUr{Nf)) 

for the zero mode component of the fermion fields 

x„W = /A*w«„(.,.). x„w = /<;S*w«t,,,,), 

where (go, ^o) are the basic elements of the zero mode subspace of the Grassmann 
algebra Q in Eq. ([T0|). As has been mentioned, due to the translation invariance of the 
vacuum field, there is a continuum of fermion zero modes, and their condensation in 
the infinite volume produces a very strong effect. Consequences of this effect in meson 
phenomenology are discussed in 0. 



5 Confinement 

Now we will comment briefly on the quark confinement produced by the field under 
consideration. Fourier transform of the two-point quark correlator defined by Eqs. ( [T0| ) 
and (ll) 

{^f{x)i!f{y))B= / dai3j2Sf{x,y) + (pert.corr.), (31) 

with S being the solution to Eq. (PT|), is an entire analytical function in the complex 
momentum plane (for explicit form of S see ^ IHl)- This means that there are no 
poles corresponding to free quarks. The other side of this peculiarity is that the Dirac 
equation for massive quarks in the presence of the background field (^ 

(iVx - m/j 4j{x) = 

has only the trivial solution = 0. Therefore, one has no appropriate field to construct 
asymptotic free states for quarks. In this sense, the quarks cannot exist as free particles 
but can propagate as virtual objects. The characteristic scale of propagation of these 
quark "virtons" is determined by the strength B of the vacuum gluon field. This 
situation can be seen as the quark confinement, for which gluon condensate plays 
the role of an order parameter. Meantime, nothing preserves these virtons to form a 
colorless composite particle by means of gluon exchange P, || . A colorless bound state 
does not feel the confining vacuum field and can be observable. A mathematically 
consistent treatment of this physical concept, especially a reasonable solution of the 
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bound state problem in terms of composite fields, requires an application of the methods 
of nonlocal quantum field theory ^, |TB|, |T7|. 

In conclusion we would like to mention the "flaws" in this picture. The problem 
about the minimum of the effective potential is not solved. Besides confined modes of 
the gluon field (in the same sense as for the quarks), free gluons appear to be allowed. 
At first sight, the gluons, longitudinal in the color space with respect to the vacuum 
field, seem to be not confined Solution of the Ua{^) problem, which is missed 

in our consideration, can come from the investigation of the local instanton-like (anti- 
) self-dual deformations of the homogeneous background field - chromons [0, rD|. In the 



meantime, these "flaws" are problems for further consideration rather than reasons to 
reject the whole physical concept. 



6 Appendix 

Here we will find the explicit form of the solution ipQ to Eq. (|1T]) corresponding to the 
zero eigenvalue: 

7mV;.^o = 0, (32) 
where the background field can be taken in the form 

= - iB^„ = ^nB^^Xu, n = ts sin C + cos C- 

Bu = B, B34 = eB, e = ±1, 
B^pBpjy B (5^1/, B^p aBp^y. 

Let us solve the eigenvalue problem for the squared Dirac operator 

(-V=^ + = ^0- (33) 

According to (^2l), the zero mode ipo has the form 

where 0o is the zero mode solution = 0) of Eq. (^31). Because of the relations ('+' - 
self-dual field, '— ' - anti-self-dual field) 

(rpuBpuP± = 0, [(JpuBf,^,'^^^] = 0, 

solution (p of Eq. (p3|) is a (right-)left-handed spinor for the ( ant i-) self- dual field, and 
it can be represented in the form 

0(a;)=^iW/(x), 75nW = ±MW, s = 1,2, (34) 
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where f{x) is a scalar function. In the Weil representation 



(1,0,0,0), 4 ^^ = (0,1,0,0) 



(0,0,1,0), 



u 



(+)t 



(0,0,0,1), 



± 5 



s=l,2 



a) 
(1) 



diag(l, 0,0,0), p!') 
diag(0, 0,1,0), Pf 



diag(0, 1,0,0), 
diag(0, 0,0,1). 



AB 



-P 



(1) 



cr P M^^) 



(35) 



The following useful relations take place 







Taking into account Eqs. (p3D-(p5|), one obtains 



(s = l,2, s' = l,2), 
/i = l,2, 

= 3, 4, s 7^ s'. 



-Vl - i-ini T e)nB f{x)=^ff{x) 



(36) 



(37) 



It should be stressed here, that for an arbitrary translation x x — xq the differential 
operator in the LHS of Eq. ( |57D transforms as [{xBxq) = x^B^j^^Xq] 



-Vl-{-iniTe)nB 



e 2 



{xBxo) 



Vl-{-iy{l^e)nB 



e 2 



[xBxo) 



which is a result of the following transformation property of the background field 



nB^yX^ nB^,{x^ - x^) = e-^^^^^^^nP^.x-^e^^^^^") - ^^^(xPxo). (38) 
This means that each eigenvalue is infinitely degenerate, since for any xq the function 



F(x,Xo) = e^("^"o)/(x-a;o) 



(39) 



is the eigenfunction with the eigenvalue ^f. 

Equation (pTf ) shows that /(x) is an eigenfunction of the operator 



n'^B 



■9 ^ 



where we have denoted 



'PVn2 



•^1 b^u 



n B 



2 B 



-'p.l/ J 
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Let us introduce the projection matrix 



then the following relations take place 



""'^■(r/ + d)Q^ir] -d) + ^^(r/ - d)Q-{r] + d). 



Using these formulas one can get 



-V^ = 2BVn^{a+Q-a + 2BVn 



V2 



Thus, we arrive at the harmonic oscillator algebra. Therefore, equation ( ^T]) and its 
eigenfunctions and eigenvalues can be written in the form 



n'^B 



2a+Q-a + 2 - (-1)^(1 t e)sign{n)\ /, = 



(40) 



ifMk,k,k, = 2Vn^B[h + h + h + h+l 



.l)«ii±ilsignH], 



aafoiv) = 0, /o = exp |-2^^} = |-^v^5a;^| 
{a^Q'a)fk^k2k3k4{r]) = {h + k2 + h + h) fk^kik-ikAv) ■ 



(41) 



The sign "±" in ( ^OD relates to the left- (+) and right-handed (— ) spinors u^^^ and 
e = ±1 for the self-dual (+) and anti-self-dual (— ) field. 

We are looking for the zero mode solutions ipo = i'~ipV ^(pQ of ( |32D corresponding to 
^^0000 = in Eq. (|40|) , i.e., we have to satisfy the condition 

(-l)^(lTe)sign(n) = 2, 

s = l,2, e = — 1, 1, sign(n) = — 1, 1. 



According to Eqs. ( 

^o(a;,Xo) 



and (|39|), (|4l|) , the zero mode solution looks as 



{n'B^fl^ 



(42) 



f{x,Xo) = exp |-iVn25(x - Xq)^ + ^nx^5^i.Xo| 



u = u 



(T) 



(e = ±1, sign(n) = 1), u = m^^^ (e = ±1, sign(n 
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The normalization constant is chosen to provide the proper normahzation of the pro- 
jector onto the zero mode subspace Vq to be calculated below. 

Relations ([36| ) are invariant under the change ± =p or s ^ s', and the zero 
mode projection operator can be constructed by the procedure being common for all 
the possibilities contained in P^). We consider the case sign(n) = 1. The projector 
operator is defined by the formula 



/^W7.V.(x)/(.,.o)Pfr(.,Xo).7.V.(.) 



47r2 

Vm (x) =d^ -iB^{x), Vu iv) =du -iB,{y). (43) 
One can check that 

j d^xof{x,Xo)f*{y,Xo) = 

V/. ix)f{x, y) = -^ (Vr}?B{x - y)f, + inB^{x - y)) /(x, y), 

i {x)fix, y)i (y) = [Q^-^ + QtA^ - y)aQlp{x - y)p] 7M^i^^7i^, 

Using Eqs. (^), we get 

'Po{x,y) = -j^f{x,y)P±, 

J d'^zVo{x,z)Vo{z,y) =Vo{x,y), 
52 

TrPo = V—. 
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